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C H A P T E R 0 N E 


INTRODUCTION 





This report continues the development of the equations of motion of a three-degree-of- 
freedpm typical section in two-dimensional, incompressible flow begun by Balakrishnan in , 
reference 1. The goal of this research is to investigate active aeroelastic control concepts 
such as flutter control and gust alleviation, although the present study is of independent 
interest as a development of unsteady aerodynamic theory. 

The design of active aeroelastic control schemes is greatly aided by the availability of 
mathematical models valid for arbitrary motions. Reference 2 discusses methods of deriving 
such models from airloads derived for oscillatory motions; references 3 and 4 present similar 
formulations for Theodorsen’s problem involving inverse Fourier transformation to obtain impulse 
response function airloads for use in convolution integral solutions of the equations of 
motion. Although these formulations are mathematically correct,, the calculations are cumber- 
some and involve functions available only in tabulated form. . Consequently, few examples of 
the exact transient response of airfoils excited by either control surface motion or gusts 
have been calculated. More common are calculations using finite state, rational function 
approximations for the unsteady aerodynamic airloads. 

In order to provide a basis for the analysis of aeroelastic control schemes, this report 
develops the exact transient response of a three -degree -of- freedom airfoil (pitch, plunge, 
and flap) forces by flap motions and gusts. The development makes extensive use of special 
time-domain functions derived from a function studied by Kussner (Sears, ref. 5). 

Chapter 2 summarizes the basic equations of incompressible, two-dimensional flow and 
Chapter 3 derives expressions for the circulation on the airfoil. Chapters 4, 5, and 6 derive 
the lift, pitching moment, and flap moment expressions, while Chapter 7 adds the section 
dynamics to give the complete equations of motion as a set of coupled, integro-differential 
equations. Chapter 8 then gives the steady state airfoil displacements, pressure distribution, 
and airloads due to flap deflection. In Chapter 9, approximations to the transient response 
model are introduced for the special cases of zero stream velocity, small time, and large time. 
In addition, a numerical solution technique is given for the solution of the general case and 
examples of the exact transient response of an airfoil at several speeds are presented. 

Chapter 10 completes the report with the development of the airloads upon the airfoil due to 
the penetration of a frozen gust field represented as a stationary Gaussian random process. 

Two appendices summarize the integral formulas and special functions used in the report. 
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C H A P T E R TWO 


THE AERODYNAMIC MODEL: BASIC EQUATIONS 




2. The Aerodynamic Model. 


We consider a typical section , of mass per lanit length m^ , extending along the X-axis 
fran -1 to +1 , with motion entirely in the X-Z plane. Let <t>(t,x,z) denote the velocity 

potential. Then 


2 2 

3 4 3 

2 r 
3x 3z 


= 0, t > 0, X < 


with the boundary conditions; 

1. Flow Tangency Condition: 



11^ Ct,x,0+) = w^(t,x) -1 < X < +1-, 0 < t 

where w^(t,x) is the downwash to be specified later, 

2. Zero Pressure-Discontinuity Condition: 

i{»(t,x,0+) - if)(t,x,0-) = 0 1 ^ X < 1 + Ut 

= O j X = 1-, (Kutta Condition) 

where Jj(t,XjZ) is the acceleration potential defined by 



4>(t ,x,z) 


9(}) 

9t 


(X» 



nie [)ertui"bation pressure is given by P = 


-P'l> • 


Following Schwarz (see Ashley [33), we seek a solution in the form: 


rl ^ .1+Ut ^ 

<l>(t,x,z) = -d/27r) I ya(t,C) Tan — — - dg -(1/2 tt) I y (t,C) Tan ^ 

X ~ 4 j w 


1+Ut 
1 


X - C 


dC 


vdiere is the circulation on the foil, and y ^"*^ 5 ^) ih® circulation in the wake. The 

Cl W 

integrals are to be taken in the Cauchy sense . 

It is easily shown that: 

Jy (t,x) -1 < X < +1 

3. ■ 

= iY„(t,x) +1 < X < 1 + Ut 

W 

and that 

= C-1) CtjXjOt) 

H Ct,x,0+) = 0, X < -1, X > 1 + Ut 


II (t,x,0-) 


d({> 

9x 


(t ,x,0+) 


f fence, in particular; 


Wt.x.ot) = j d? -1 < X < +1 

= > dC + j ^ Y^(t,C) dc), 1 < X < 1 + Ut 


Again, even though we cannot differentiate with respect to time inside the integrals, we have 


to 



by considering corresponding differentials. 


3 

3 t 


(J>(t,x,0+) 


3 

3t 


4)(t,x,0-) 


Hence the Boundary Condition 2 yields: 


Uy^ (t,l+) + r»(t) = 0 

And Kutta Condition yields: 

Uy^ Ct,l-) + r'(t) = 0 


Hence, 




Mcreover, we have, from boundary condition 2 


In particular 


and is finite. 


Hence 


Putting 


or. 


or. 


Hence 


" (t,x) must satisfy: 
w 


Y^(t,x) = F(t - J ) 


X = 1 + , we have 


= «t-i) = 


f'(t) = U F(t - 1/U) 


F(t) 


r'(t + 1/U) 

u 


Yw^t,x) 



+ 



-r'(t) 

u 


t S 0 





Now 


w (t,x) = 

a. 


1 Y^(t,y)dy ^ ^1+Ut y,^('t,y)dy 


2ir J_i X - y 2 tt I . 


w 

X - y 


Substituting for y (t,x) , we get; 

w 


1 Y^(t,y)dy 


- f 

2 tt 


X - y 




1+Ut r v(t + )dy 


^ 

X - y 


= w(t,x).^ f 

a 2 ir J - X - 1 - 


r ’ (t - g)dg , . 

W~ V -1 < X < 1 


Now, this integral eq\aation has a vinique solution, since y^Ct,!-) is required to be finite; and 
given by: 


xs 


(t,x) = ^ I 

It/ 1+x I 7t-y X - y ^ 2ir J_^/ 1-y x-y J ^ y - 1 - Ug < 


This is our basic integral equation. 
Let 


.1 /rr- 

H(g,x) = I = . . dv s (_Tr) 1 /z+1 , 

J_i / 1-y y - 1 - Ug x-y y ^ (x - z) / ’ 


z = 1 + IJg 


Tlien the equation beccmea; 




1 ^ 

l--y X - y 



giving the total circulaticn on the airfoil as 


fl /tT 1 /*t 

r(t) = 2 J \/ w^(t,y)dy - ^ jf H(a) r'(t-a) 


- 1 /* /£i 

■n J V 1+x 


1-x 

H(o,x)dx 


(see Appendix 1.) 




CHAPTER THREE 
CIRCULATION ON THE FOIL 




3. 


The Circulation on the Foil. 

To proceed furtlier, we endow the foil with three degrees of freedom. With reference to 
Figure 1, let h(t) denote the plunge coordinate, a(t) the angle of attack, and 3(t) the flap 
deflection. Then the downwash w (t,x) is given by 

3 . 


where 


w^(t,x) = Ip 2 ^(t,x) + U 2 ^(t,x) 


3x "'a' 


z (t,x) = 'T\(t) - Cx - a) aCt) , 

3 


-1 < X < c 


Then 


- -h(t) - (x - a) a(t) - (x - c) 3(t) , c < x < +1 


w (t,x) = - h*(t) - (x-a)a'(t) ~ Ux(t), - 1< x < c 

a * 


= - h'(t) - (x-a)<^’(t) - Ua(t) - (x-c)3‘(t) - U3(t), c < x < 1 



g Hence we have for the circulation: 


r(t) = { y it ^x)dx 

/-I "" 

2 r- /!-?. dx 

'-1 ,/I+^ 


• ^ 

^-1 ^ 


w^(t.c)dc -5j H(a) 


(- 2 ) 


C w (t,a)d0 

^-1 ^ 




\ H(a) r*(t-a)da 
•'o 


Hie last line follows since 




-1 vl+ic 


2 

7f 




/I7 


-Ddx = - 2. 


Also 


£. /l-x dx r^‘ 

'-1 A+x '-1 


w_(t,i;)dr, 

A -T ^ 


f' 


A+F 


' (t-a)da 



= 2 


n 1 

(h«(t) + ( 7 ;-a)a'(t) + Ua(t))dC + 2 [ 


/UC 


((C-C)fi'(t) - 


= 2tt h' (t) + 2tt Ua(t) + (fr- 27 ra)a 


' — 7 —1 2 ^ 

'(t)+ (2^-0^^ + (l-2c)cx)S c - c/l-c )3*(t) 


+ (2gos“^c + 2 / 1 -c )U 3 (t), 


Hence 


r(t) = 2 irUa(t) + ( 2 cos'‘^c + 2 » 4 -c^)Ui 3 (t) 


+ 2 ir h’(t) + (Tr-2TTa)a‘(t) + ( 2 ^C? + (l- 2 c)cos“^c - cvO-c ) 3 ’(t) 


f (i_ J^)r«(t-a)da. 

•^0 ^ 


Let x(t) = 


h(t) 

aCt) 

3 (t) 


I- U 3 (t))d^ 




Then we can write 


r(t) = CB,Z(t)] + 


where 


B 


But 


r(0) = [B,Z(0)] 


f (1 _ )r»(t-a)da 
Jo »CT 


0 

2irU 

(2 cos~^c + 2v4-c^)U 
2TT 

TT - 2ira 

(1-2c)cos'"^c + (2 - c) 





fjence 


0 = [B,Z(t)-Z(0)]- f r'(t-a)da 

■’o 


Now for Res > 0; 


^-so /z-l-'l 


da = 


= ^ S. 


00 _ i */z+r 

U dz 


»^z-l 


Now 


Let 


+ s ^ ^ 


1 U f 

= tr - 


/z-1 


sz 

u /z+r _ 

e dz = 

.CT 


sz 

U (z+1) 




dz 


KflCs) 


CO 

1. 


-St dt 


/t^-] 



K> 

JSJ 


Hence 


-St t dt 


1 






SZ 

\F 




Here Kq(s) is the modified Bessel function of order zero. 
Hence 



•^da = (Kq(|) - Kg(s/U)) 


✓z^T 


U 


and hence 



e" ^‘^r*(a)da 


UCB,L(Z;s). - 2v0ljg-s/U 

’ s 

(Kq(^) - lys/u)) 



v;hexe 


L(Z;s) = ( e ®^'z(t)dt 

^0 


inverse Laplace transform of 

^ 1 1 

^ KflCg) - k;(s/U) " r ^-sa ^ 

h vCT 

by 

c^(t). 

‘Ilien 

I‘‘(t) = 1 CT(t-a)CB,Z(0)]da 

^0 ^ 


Let us denote the 


U e 


-s/U 


ro 

CO 



K) 

■P 


and 


r(t) = j c^(t-a)CB,Z(a)]da - C2(t)[B,Z(0)3 + [B,Z(0)] 


wheire 


C 2 (t) = I c^(s)ds 


Finally we note the series expansion due to Kussner-Sears [53 for c^(t) 


Ut < 5 


G^(t) = U ^ i (Ut) 


- 1/2 


^ ^ (Ut)3/2 




26380 


Ut > 5 


U 


c (t) = { 1 --^ 

(Ut-1) (Ut-1) 


^ 4Log C2Ut - 2) 
^ ■■“■(iJt-1) — ” 


54 Log (2Ut - 


2) (14 + i - 3 tt^) 

■' ZrT 


(Ut-1)' 


(Ut-1)' 


+ 18 (Log (2ut - 2))' 
(Ut-1)^ 



c^(t;U) = U c^(Ut;l) 
c^diU) = c^CUt;!) 


Figures 2 and 3 present plots of c^(t) and C 2 (t) calculated using the above approxinations . 


ro 

cn 




CHAP TER FOUR 


CALCULATION OF THE LIFT 




4 . 


Calculation of the Lift P 


Vte have 


P = I P(x,t)dx 


= (-P) [ + 


Y^(t,y)dy dx] 


= (-P) Cu ret) - 


al’ { xY^(t,x)dx + r*(t)3 


Now 


Now 


_3 

\ 3CY 

(t,x)dx = 1 

.1 2 

i 

■^X 

dx ( 


dr a 

St 

J_1 


Ll^ 

/l+x 

J 

-1 /1R 

X— C 



a 

“ at 

iO 

1 

X 

2tt 

H(a,x)dx 

r'(t-o)do=T^ + T 2 


/Tk 

jIw (t,C) 

- -f. '■ dC ^ 

= 

/Ik 

1 

(-h“(t) - (C-a)a"(t) - U 

J-1 


^-4 


'-1 

/R 

x-C 


to 


i: 


vtk 1 


(- (c-c)3”(t) 


•(t))dc 


U 3'(t))dc 



03 

O 


= +7T h”(t) - (air - IT - irx) a'*(t) + Utr a’(t). + (cos~^c + 0”(t) 

-1 / i~^x 

+ (x-c) (cos c + / Log ||) 2"(t) 

+ U 3'(t) (cos ^c + — ■ Log j I ) 

VT=x 

where 


Log hi =iog Ai-c)(i->x) ^ /g-tcia-x) 

|Al-c)(l+x) - Al+o)(l-x) t 


I? ‘-g N 




/L-x 


1 

x-c 


and vdiere we note that 


1 * — X 

( x(l+x) - dx = 0; 

f-1 


/1+x 


J ^ X log II dx - i c J\-o \ I x(x-c) log I* I dx = i (1-c^)^^^ 




x(x-c) 


!^^dx = I(H-c); ( 

^ 2 y_3_ 


,, ^ if 

X dx = - ^ 



hence 


= - IF h"(t) 


+ (air) a”(t) - U ira'(t) + (c cos~^c - )^) 3 ”(t) + U(c 


A-c? - oos“^ 


Next 


T = 
2 


t 1 


a [ r r 2 ^ ^ . (_ 1 ) H (a.x)dx r'(t-a)do] 

^0 3_i ^ 2ir 


= [ J (1 + Uo - /(l+Ua)^ -1 )r*(t-o)da] 


= + r*(t) + 


at 


( (Ua - vW + 2 Uo) r'(t-a)da 

'O 


Hence 


P = (^p) {Ur(t) + Ujra‘(t)- U(cX-c^ - cos"-^c) 3 ' (t) +ir h‘'(t) - aira ”(t) 


-1 


- ( I + c cos"^c - 3"(t) - ^ CgCt-a) CB,Z(a)3da} 


c) 3 * (t) 



CO 

K> 


where 



The function c^Ct) may be calculated numerically using the approximations given in Chapter 3. 
Figure 4 presents a plot of CgCt). 



CHAPTER FIVE 


CALCULATION OF THE PITCHING MOMENT M^ 




Calculation of the Pitching Manent M 


We have: 


X 

= (-p) { ^ (x-a)U Y^(t>x)dx + 


(x-a)dx 


Y^<t,y)dy) 


3 f ^ (x-a)^ 


(1-a)^ 


{ J ^U(x-a)Y^(t,x)dx - ^ J ^ Y3<t,x)dx + i±J 


r*(t)} 


I (x-a)Y^(t,x)dx = - a F(t) + j xy (t,x)dx 
J _1 a ' -1 ^ 


( X Y (t,x)dx = ( f • X • { f 

^-1 J -1 "" /IT^ J -1 


/Ik 

/Tk 


w (t,C) 


1 

- 7j— I H(o,x) r'(t-o)do} dx 
2ir Jn 


j 2 /^ ^ ^ f dc + ( (d+Lto) - */(jV+2UQr)r'(t-o)da 

»-l ^ ✓I+x 1-1 /Tk Jo 


. t 

= + (r(t) - r(o)) + J Cgd-o) [B,z( 0 )]do 


CO 

cn 



. - 2 2 ^ ^ [. h'(t) - (;-a)a(t) - Ua(t)l • ^ 

1 ■ U /m ^ -1 ^ 


+ i5li [ - (c-c)e’(t) - u 3(t)] ^ } 

3 c 


^ 2^ X • {Tr(h*(t) + U a(t)) - tr(a-l-x)a’(t) + U (cos ^c 

-1 ^ ^ 

+ C(cos"^c + A-c^) + (cos”^c + log I I )(x-c)] 3 ’(t)} dx 

/E=x 


= - IT (h*(t) + UaCt)) + TTa a’ (t) + U(-cos"^c + c/C?) 3 (t) + Cc cos'^c - | (2 + c^)/C? 33 '(t) 

Next 5 the **non-cii?culatory” terms (that is, terms not containing F(t) ) in, 

si \\ Y^<t,x)dx 

i <x-a)^ In (h"(t) + U o'(t)) - x(a-l-x)o"(t) + U (oos'^c + log 

)-l ” 

+ [(cos'^c + A-o^) + (cos"^o + log |•|Kx-o)B"(t)} 


2 

= (1 + + a)TT(h"(t)+Ua'(t)) + tt(|- - y - )a"(t) + UC(y + a + a^)cos"^c 


+ (i + 2 a + ^ - 2 ac) 
6 d 


X-c^33' (t) [(cos + X-c^)(i + a^ + a) 


+ }- (4ac^ - 6a^c - - Ma - |)- cos"^c(a + f~ + | J ^ * ac)33"(t) 


12 


Hence 
M 


a = (-p) {-aUT(t) + U [F(t) - r(0)3 + U f c (t-a) [B,Z(a)3da 

JQ ^ 


+ U[- irh'Ct) - Ud(t)ir + iraa' (t) + 


U (- cos ^c + cJ\- c^)3(t) + (c GOs”^c - i (2+c^)X-c^ 


+ r'(t) - (i + a^ + a)7T (h"(t) + Ux'(t)) - ir(i - | - a^ - |-) a"(t) 


- U [(^ + a + a ) cos ^c + (^ + 2a^ + — 2ac)/^ l-c^33*(t) 


CaJ 


- [(cos ^c(^ + ^ ^ - ac) + 


X-c^ (i + 


^ + a. + ^ a c — 


+ 4ri ( [i_ f (x-a)^ H(a,x)dxjr’ (t-a)da> 

*2 1 ^ ^2 ^ 


)3'(t)3 



where the factor in squcire brackets in the integrand in the last term is 


[ (x-a)^ H(g,x)dx 

TT J -1 


= (i - 2 ^ + 2 + 2a (2 - - i) 

It is also possible to split llie circulatory and non-circulatory terms in yet another way. Ihus by 
not splitting 

J ^ (x-a) Y^(t,x)dx = - aT(t) + j XY^(t,x)dx 


but directly calculating the left side, we have: 


J ^ (x-a)y^(t,x)dx = 


j I (x-a) {Trh'(t) + irUaCt) + Trd-x-a)a'(t) + U(cos"^c + log | j )3(t) 

-1 /1+x /[lx 


1 + C(oos ^c + /C?) + (cos"^c + log j |) (x-c)33'(t)} ^ 

/Px 


= - (l+2a)TfUa(l) - (l+2a)nh' (t) + 2a^Tra'(t) + u ((-l-2a) cos”^c + (c-2a)X~c^)3(t) 


+ C(cos“^c + (-l- 2 a) + 2 coa"^c (ac + ~ + i) + | /C? ( 1 -c^ + 


3ac)'J3’(t) 


T 4 = + 2 tf 


( tz - »4^-l + a - a] r*(t- 0 )da 


This yields; 

= (-p) {U^(-l-2a)Tra(t) + U^C(-l-2a)cos“^c + (c-2a)X-c^]B(t) 


- U(l+2a)nh'(t) + 2Ua^ira'(t) + UCj Jl-c? (1-c^ + 3ac) + 2cos~^c(ac ^ ^ I*) ~ (l+2a)(oos~^c + /{-c^)]3‘(t) 

- (i + a^ + a)n U a'(t) - U [(^ + a + a^)cos"^c + 4 + 2a^ + ^ - 2ac)X-c^ 33* (t) 


- (| + a^ + a)ir h"(t) - Tr(| - - |-)a''(t) 


- [(cos"^c + + a) + /C? (4ac^ - 6a^c - - 4a - J)- icos~\Usi + ^ + | + | + aVac)33"(t) 


(1-a)^ 


r'(t) + U ( (z - A'^-1 + a - a)r*(t-a)da 


I ^ J [j - z^ + zJ?^ + 2a (z - 


+ a^ - l)]r'(t-a)do 

.CT 


CaJ 

CO 



(-P) {U^(-l-2a)TT a(t) + U^C(-l-2a)cos"^c + (c-2a) A-c^]6(t) 


‘-U(l+2a)Tr h'(t) + Lhr(a^ - a - i)a'(t) 


2 2c^ 


+ U[(- ~ + 3ac - 2a - 2a'^ - ^ )A-c^ + (2ac + 


(2ac + c - j - 2a - a^) cos ^c]3'(t) 


(| + a^ + a)ir h»(t) - | - ^)7ra»(t) 


- - - ac)cos“^c + A-c^ (4 + + I- a + |- ac^ - - 


I c. 

12 ■ 24 


)]0” 


2 

2 "" r'(t) + U f (z - /z -1 + a - a)r'(t-o)da 

■•o 


+ ±. .2 ^ /2~ ^ 
^ 7?r I T - z + z/z -1 + 


The last two lines may be written as 


2a(z - /z^-1) + _ l)]r*(t-a)da } 

/z-1' 


r»(t) 


Itl f f*- 

2 — + U 1 Cjj(t-o) [B,Z(a)]da + U a (l-c„(t-a)) [B, (a)]do 

“^0 n ^ 


dt 


1 

t 4 j (l-C 2 (t-a))[B,Z(o)]da + ^ r(t) 


1 ^ 2 
^ 2 G^(t-a)[B,Z(a)]da + a f c (t- 0 )[B,Z(a)]da + [B,Z(t)]} 

4 0 , -»0 


[1 c (t-a) + U Cj^(t-a) + Ua(l-G 2 (t-ai)]CB,Z(a)]da 
Jo ^ 

+ 1_ [(i- + ^) - i c„(t-a) + a cAt-a) + ^ c^(t-o)][B,Z(a)]da 

dt Jq 4 2 *4 2 O 

Where 


c^(t) 


= J c^^Ct - a) (1 - Ua - ^ + 2Ua^ 


da 


= C2(t) + CgCt) 


Cg(t) = 


J Cj^Ct - o) ^(1 + Uo) yj * 2\Ja - (1 + Ua)V 


da 


Plots of Cj^(t) and Cg(t) are given in Figures 5 and 6, 




CHAPTER SIX 


CALCULATION OF THE FLAP MOMENT M, 




Calculation of the Flap Moment M.. 

P 


We liave: 


4. 

Mp = ^ (x-c)P(x,t)dx 


- (-p) { J U(x-o)Y^(t,x)dx <■ jI- J (x-o) dx J Y^Ct.y)!^} 

= (-P) ( 5 U(x-o)Y^(t,x)dx + g| Y^(t,y)dy- 5 | i2=|l' 

C -1 n 


^ (1-c)^ 


= (-p) { U J (x-c)Y^(t,x)dx + F'(t) - 


1 

2 at 


9 

\ (x-c)“ Y^('t,x)dx} 


Now 


*-i 


w (t,C) 

ci 


dC 




-1 


✓T+T (-h*(t) -• (C-a)a'(t) - Ubt(t))- ,^ /I+^ (-(c-c)8'(t) 

3 “ — 


'c 


(x-Cl 


Y^(t,x)dx) 

a 


U8(t))dC 



-p 

CD 


Tl'ie term 



F^(c) 


dC = F2 (c) 



3 




, vCT: r -1 k 5* , \ -1 ^ 2(z-i)(z-c) -1 /rra*^+T 

= (- 2 ) Coos c - / 1 -G - (z-c) cos c + Tan — 




v4^ 




Hence 


^ (x-c)Y^(t,x)dx 


/L-x 


j ^ (x-c) (frh’Ct) f ir(x-l)a' (t) ■{• aira' (t) + U»ia(t))dx 




- (F 2 <c) - c F^(c))3'(t) - U F^(c)3(t) - ~ J H^(a)r'(t-a)dcf 


= 2(h*(t) - a a'(t) + U a(t)) ((1 i- J\-o - (c ^-) cos''^c) 


2 

+ C(|+ |-) A-c^ - oos"^c]a'(t) 

- (F 2 (c)-g F^(c))3‘(t) - U F^(c)3(t) 
1 ^ 

"H Hl(a)r'(t-a)d0 


-F 

00 


Next 


5 (x-c)\^(t,x)dx 

c 

.2 

yl+x 


— C (x-c)' 
TT J _ 


(irh'Ct) - aira'Ct) + U7Ta(t))dx 


i- 2 C (x-c)^ (xfl)dx a’(t) 


- (gj(o) - ogj^(o))C’(t) - U gj^(o)6(t) - |j [ H2(c)r'(t-<i)do 


where 


H2(o) 


glVC) - - 


1 

( (x-c)' 

J. 


c 


So 





l-x 

1+X 

H(a, 

x)dx 

-/l-x 

dX 

1 

( 

A+i;' 

■/Id^ 





dx 


/I+? 



c 



x-c 


ciC 


x-C 


Plots of H^(t) and H 2 (t) are given in Figures 7 and 8. 
Now using formulas (23-24) from appendix 1, we have; 


■4 


First Term = 2 hj^(c) (h'(t) - aa*(t) + Ua(t)) 


Second Term = + 2 h 2 (c)a'(t) 


Next 


a 


,(o) = I 5 (x-o) 


2 


A+x 


H(o,x)dx 


= (-2) ( \ 


^ (x-c)^ A-x A^•l 




/H’X 




(.2) ^ + x-2ctz]dx 

/z-1 ' C Atx X z 


(.2) 2^ [(z-o)2(-cos-lc t 2 Tan"^ ^ 




/T7’ /D+HTAa 

r Z -1 


+ (z-2c)(cOS ~C - A~C^) + (*4~C^) - ~ COs’”^C - 


-p 

to 


(-4)(z-c)^ Tan~^ - 2 C(z-2c-(z-c)^ - i )cos‘*^c + (l+2c - 

/ttc)'v''za: vCT 


N>jO 



Hence putting it altogether we have; 


■TFT 


U((2+g)X-g^ - (2c+l)cos ^c) (h’(t) - a a'(t) + U a(t)) 


UCCj + - 


c cos ^cDa’Ct) 


- U (F 2 (c) - c F^(c))p'(t) - Fj^(c)0(t) 


2 

" H^(a)r'(t-a)da + r»(t) - h (c) (h"(t) - aa'*(t) + Ubt‘(t)) - h.(c)a''(t) 

0 ^ 


^ +1 g^(c) 3 '(t)+i (i^) ^ J%2(a)r'(t-a)d0 

= (-P) {u 2 c( 2 +c)/I? - ( 2 c+l)cos~^c]a(t) - f^(c) 3 (t) + U [( 2 +c)^I 7 - ( 2 c+l)oos"^c]h' 


(t) 


2 ( ) 

+ U C(| t c t - a+3o)oos"^o - hj^<c)]a'(t) + U - (fj(o) - ofj^(o)]6'(t) 

- h,(o)h”(t) * (ah^(c) - h2(c))a”(t) . | (g^CO - o g^(c))S"(t) - H^(„)r.(t-a)da 


* I 5 ? it 1„ HjCOr'tt-Oda ) 



The functions f2(c), gj^(c), g2(c) are calculated in Appendix II. 

h, (c) = f (x-c)^dx 

■'c /L+x 


,li2c+2c^v -1 1 /o 2 ^ , I 2 

= ( 2‘- -) cos c - g- ( 2 c + 9 c + i|)A-c 


lu(c) = f (x+lXx-c)^dx 

•'c A+x 

= <■■ Mo) 


oi 




CHAPTER SEVEN 


EQUATIONS OF MOTION 





7. Equations of Motion. 


h(t) 

Let x(t) = aCt) 

e(t) 

Then the typical section equations of motion with trailing-edge control can be written: 

M *x + B X + K X ~ + Gu 
a s s m 

s 

wher \2 the subscript s stands for structure, and 

'a ^3 

, Cr^+xg(c-a)] 

Cr>^+Xg(c-a)] 




cn 

cn 






p = p + p 

c nc 


M = M + M 
a a,c a,nc 



where the subscript c stands for 'circulatory' and nc for non-circulatory. Now we can write 



wjiere the subscript 'a‘ stands for "aerodynamic" 


oi 

CO 


0 


0 


0 


K = [lil] 
a m 
s 


0 U^(-.l~2a)ir 


irU 


% 

C(-l-2a)oos~^c + (c-2a)«4!^3 


[o U^C(2+c)X? - (l+2c)oos“^c] 


- U‘f^(c) 


Note the difference in the 2nd row corresponding to the pitching ironent terms. 


B 


= C^] 

m 


0 


-irU(l+2a) 


Ut(2+c) 




(H2c)cos~^c] 


irU 


7fU(a^-a- 


UCc 


Xj. 


oob”^c] 


U {(- ^ + 3ac - 2a 




UC(|- + c + ^)vC? 

- (l+3c)cos“^c - hj^(o)] 


1 2 1 
- ( j + a + 2a - 2ac - c)cos c 

^^5^ “ <fo(c) - c f;(c»] 

O * X 


} 



M 

a 


- 


C COS 


“^c - ~ ^ 




(i>£ 

'“8 ^ r 


“I-ac 


2/i_c 2 (1 + a 


+ |:+ ^ 
3 6 


a h^(c) - h2(c) 


^ (§2(0) - c gj^(c)) 


Next let; 


Z(t) = Col. ( xCt), x(t) ). 


Recall that: 


2TTU ^ 

(2cos"^c + 2 J (l-c^))U 


(l-2a)ir 


(l-2c)cos*c +(2-c) 4 (1-c^) 



o> 

o 


Finally ; 



U ( c (t-o)[B,2.(0)3da + U CB,z(0)] - — | c (t-a)[B,Z(a)]da 

In 2 at Jj. J 


[—3 f c (t-o) + U c, (t-c) + U a(l-e„(t-a))3[B,Z(0)3da 

IDg )0 >4 1 t 2 

,t 2 , 


^ it \ *^^5“ ^ ? c^Ct-a) + a CgCt-a) + | Cg(t-a)3[B,tCo)3da 


~ I Cg(t-a)[B,'Z(a)]da+ j J c^(t-a)CB,i(a)3da 



M„(t-o)Z(a)da 
0 ^ 



Mo(t-a)Z(a)dcf 
0 ^ 


+ U[B,Z(0)3 

0 
0 



>wl ,e 


OjCt) 


I Cj^(t-o) (Ife - Aw * U‘V )da 


= j c^(t-a) (1 + Ua ~ Ai^a^ + 2Ucr)da 


= + Q^(t) 


5(t) = [ G^Ct-a) (<1 + Ua) Aj^a^ + 2Ua - (1 + Ua)^)da 

-»n -*• 


r"C 

= j H^(a)da 


c? (t) = f c,’(t-a) H.(a)da 
•’o ^ 


Plots of Cg(t) and c^(t) are given in Figures 9 and 10. 



Oi 

to 


The natrices M^Ct) and M^Ct) are 



A 

Uc^dOB** 

|^|c^(t) + uc^(t) + au(i-C 2 (t))j b" 


-u 

2ir 


Cg(t) B** 



-C 3 <t) B 

2 

1(5- *k)-k * K'^>] 


1 

•fir 


c^(t) B 


The equations of motion can thus be written finally as: 

t 


MX + Bx + Kx = Gu + Mx + Bx + Kx + 
s s s a a a 


(■ M, 
•'0 


^(t-a)Z(0)da + — 
^ at 


T. 

j M3(t-o):^(a)da 


uB*Z(0) 

0 

0 


( 6 . 1 ) 




Or, 


Z(t) = AZ(t) + Hu(t) + [B,Z(0)3 [ V ] + 

V . 1 

I (M -M ) ^ M^(t-o) Z(a)d0 
Jo ® ^ ^ 

+ 

t (M -M ) ^ 
dt s a 

t 

1 M3(t-a)Z(a)do 




- 



Where 


o 

I 

3x3 

-(M^-M )~^ (K -K ) 
s a s a 

-(M -M ) 
s a 


-1 


s a 


H = 


o 


(M -M ) -^G 
s a. 


V = (M -M )"^ 

s a 




a> 

CO 


Circulation: T(t) = e^(t-a) [B,Z(a)]da + [B,Z(0)] 


'0 


Since ^ Z(0) = CB,Z(0)] 


j M^Ct-a) Z(a) do = M^C") Z(t) - M («>) Z(0) + f M ( t-<j) Z 
0 ^ i 2 


Z(a) da 


we have finally: 


B ) Z(t) + 


M 2 (t-a)Z(a) do 


(M f M 

So. I 

0 


+ H u(t) 


2(t-a)Z(a) da 


( 6 . 2 ) 


UCc^Ct)-!) B 


M^Ct) = 


(q- (c^(t) + U c (t) + a U(l-C_(t)) 


(?) 


- !Lcc(t) B 





CHAPTER EIGHT 


STEADY STATE RESPONSE TO A STEP FLAP INPUT 





8. STEADY STATE RESPONSE TO STEP FLAP INPUT 

In this section we shall calculate the steady state response 
to a step input deflection of the flap? assuming system st ab i l ity. 


8 .1 Steady State Structural Response. 

We begin wi1di the calculation of the steady state 
structural response. The dynamic equations are: 


M V + B X +K X = Gu * M_ *x + B_ x + K x + M2(“) Z(t) 
s S S a* 3. a ^ 

t • 

+/M2(t-cr) Z(a)da 
o 

-— • / MgCt-cr) Z(a) dcr 
o 


M2(t) = M2(t) - M2(“) 

Hence the steady state response obtained by setting all time derivatives to 
zero and initial conditions to be zero, readily becomest 

(Kg . K^) x(«) U B* Z(«) = G u(«) 

In particular , taking u(<» ) = 1, we have: 


69 





♦h 


(- ^) U c^(«) (2irUa^ + (2 cos“^c + 2 A-cM ) 
s . ^ 


.-1 


U^.(cos“^c - c vC?)]3 


2 > 


^ ^"’ *=2^**^ C(-a)(2TTUcc^ + (2 cos'’^c + 2 v4~c^)U3 )3 


2 2 2 2 
Or*, finally: 


(-P) 


m 


CU" ((2+c)/l-? - (l+2c)cos"'^c)a - U^f 


U C2(«) (27rU)a^- (^) U^C2(«>)(2oos“^c + 
s s 

(|^)ula + a (^) C^(~)27ru2]a„ 


2y4-c^)3 


= 0 


' -2 


+ C(+a) (^) U^C 2 («) (2cx3s”^c + 2/C?) 
s , 

+ C^)U^cos“^c - ci/^-c^)] 3^ 

ill_ CO 


= 0 


^ ^ ^ ^nT~^ ^ C(2fc)v^-c^ - (l+2c)cos ^c)3ct 

s “ 

+ {[^3 [-u2fj^(c)3 ♦ r-2o.2}8„ 


-22 
■ 3^3 


(03^] 
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Letting 


(;—) • 2ir 


r (jj - U^(2aTr+ir) 

a a m - 

s 


(^) C(2+c)XI? - (l+2c)cos"^c] 


C^)U^ (2cos~^c + 2 A-q^) 


C;~)U^C(c-2a)A-c^ - Cl+2a)cos~^c] 


§- C2i/C? cos~^c + Cl-G^) 

^ 9 . 

- (l+2c) (cos 


2 2 

* % 


We have, for the steady-state values in general. 


X («) = G u(<») 
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This yields for the data of Table 1 


< 2 . 5 ) 10 ^ 0 0 


"V 


"0 

0 ( 2 . 07 ) 10 ^ ( 6 . 25 ) 10 ^ 




0 

0 ( 2 . 6)10 ( 6 . 11 ) 10 ^ 




. 5.625 X 10 ^ 


Or 



- 

- ( 4 . 0 ) 10 ^ 

(- 7 . 9)10 ■ 

(- 6 . 9)10 


"0 


■- .39 “ 

00 

Z 

c 

( 4 . 87 ) 10 "^ 

(- 4 . 98 ) 10 “^ 


0 


- . 28 ' 

8 

00 “ 

r 

^ 0 

(- 2 . 12 ) 10 “^ 

( 1 . 65 ) 10 ”^ 


.( 5 . 625 ) 10 ^, 

1 

_ 0 . 93 _ 
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steady State Pressure, Lift and Moments 

We can also calculate the steady state pressure distribution due to 
flap deflection directly from the aerodynamic equations. We have 

Steady State Dcwnwash 


w^(*,x) = 

ol 


- Ua, - 1 < X < c 


= - Ua - U6^, c < X < 1 

OQ ao ^ 


steady State Circulation on the Foil 


Ya(“,x) = I 4^ 

® '' A+x 




A+C W (oo,;;) 

Ol 

A-? X-? 


d? 


(2Ua) - (2UB) - ^^3:^ 


A+? . 1 


/l+x 


TT A+X 


A-? 


X-? 


dj;- 


^ (2Ua) . iMI { ccs-^c - 

A+x A+x A-x 


)< ^(1~g)/1+x + /(1+c)* /1-x 
/(lie) A+x - /(1+c) /1-x 


A^ + (2US) ^g-1^ ^ (2US) 

'7T I 


/l+x 


(2Ua) 


/l+x 


ir 


Log 


/(1-c) /l+x + /(i+c) /1-x 


/(l-c) /l+x - /(l+c) /1-x 


Note that Y (*jx) is discontinuous at x = c. Note also that 


Y («, 1) = 0 = r’(«) 

3, 

Figure ll is a plot of the circulation, Y _(”5 x) , on the airfoil due to a unit 


step change in flap deflection. 
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Steady State Prassure 


?(»,x) 5 (-p) U Y (»,x) 


Steady State r.ift 


L = 




P(«,x)dx 


s (-p) U 


f Ya^*j3c)dx 

J .1 a 


C-p) U^{ -L 28 


20 cos c 


+ ^ f’ 


x4=s- -i= dx } 


x-c 




= (-fi) i/ (2m + 28 oos'^o + (28) >d-c^ } 


Steady State I'fcrcerrt M 


a 


~ ^ J Y^C"jX)dx 


= (-p) l/ {[(-IT -: 


2arr)a - ((l+2a)cos"”^Gj8 + ~ 

TT 


1. 


(x-a) A 


2 A 2 


-<2 . 1 


x-c 


/mmmmmmm 


dx} 


= (-p)l/ {-<2a+l)ira + 8 C-(2a+l)cDs"^c -t* Cc-2a) A-o^ 


]} 


74 


Steady State Mbment Mg 
-1 

Mg = (-p)U I (x-c) Y («»,x)dx 
= (-P) C(2+c) (2c+l) cos”^c]a 

+ C ^ (2+c) X-c^ - ^ (2c+l)cos ^c] (a3s”^c)]8 

= (-p)U^{ C(c+2)A-c^ - (l+2c) cos"^c]a 

+ (i) CGS’^C - (l+2c) (cos"^c)^ + (1-G^)I6} 


Finally noting that the steacfy state response to step input (xaiit flap 
deflection y is determined by: 


K_x_ 
3 00 




+ Gu 


we obtain 


1 

w 

0 

- 

0 


p 


0 

h 




■^S 



0 

2 2 
r (jj 

0 

X = 

M 

a 


0 


a a 


00 

"^s 

+ 


0 

0 

2 2 


i 


2 2 
r«£»)« 

- 



s 


B B 


leading to the same equations as before. 


In particular, for the data of Table 1, we have: 


” p ” 




2 






“h • 


- 975,00 








a 

S 

kx - Qu 
s ^ 


r^oj^a 
a a « 


- 700 

, 





2 2.- 



_“s_ 






- 38.25 


A plot of the stsatfy state circulation on the foil is given in figure 2. 
Nora the logarithmic discontinuity at the hinge point x=c. Note that 
the Kutta condition is clearly satisfied at the trailing edge. 
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9. CALCULATION OF TRANSIENT RESPONSE. 


In this section we indicate our technique for* calculating the tfansient: response of the 
three-degi*ee-of-freedora typical section to a step input to the flap. 

For this purpose it is convenient to rewrite the basic dynamic equations in the 'integrated* 
form, and taking Z(0+) to be zero, without much loss in generality. 



Integrating the equations of motion: 


Z(t) = f AZ(a) 
0 



o 




da + 

f (M -M ) M„ (t-a)Z (a) da 
J Q s a 2 

+ 

(M -M )"^ 

S 3 , 

f M,(t-a)Z(a) da 
■'o ^ 


+ Ht 


We rewrite this equation in terms of Z(t): 


. A 

I Z(a)da - / AZ(a)(t-a)da - 

Jn Jq 


o 

ft ^ 

J M^(t-a)Z(cf)da 


o 


-1 

J M^(t-a)Z(a)da 


vdiere 


M.. 


ft 

,(t) = I M_(a)do. 
*' 0 ^ 


Hence, vie need to solve , the convolution integral equation for Z(t): 


f M(t-a)y(a)da = Ht 
•' 0 


; (y(a) = Z(a)) 




00 

ro 


M (t) 
2 


MgCt) 



UC 2 (t)B 




h- c^(t) +ucj^(t) + au(l-c. 


B* 




-C 2 (t)B 


■(-t) 


|^+ . 1 

\2 4 y 4 


c^Ct) + aCgCt) Cg(t) |B' 


5?=7‘«B 


■St 


c^(t)dt = 


'/ 


g J g-Sff /2+Ug 

0 /ua 


da 


1 


3 


,4-% - . 1 t 4-A-4 1 r 4-^4 4- 5 , . x3/2 


, ut < 5 


c (t) = f c,da 
^ 0 


j(t) - J" Cj^("t— o) (c»Cf — + 2ua)da 


Gj^(t) = C2(t) + CgCt) 


CsCt) = t /„ c^(t-a) [(1 + u a) +2uo - (l+ua)^ ]da 


Cg(t) = J c^(t-o) I^(a)da 


c„(t) = f c (t-a) H-(a)da 
/ j 0 ^ 


00 

00 





OlCM 


0 


Solution for zero stream velocity : u = 


It is of interest to exandne the limiting case when U=0. We have: 


M2(t; U=0) = 0 


M3(t;U=0) 


(-P /ffig) 


((a^/2 + hi) B* 

(1/4tt ) [(l+2c^) cos”^c 


- 3c B* 


(Let us denote the matrix cai the right by q) 


K = 0 

a 


B =0. 
a 


Thus the dynamic equations for the case when U=0 can be written: 

0 


Z(t) = A Z(t) + H u(t) + 


s ^ 3^0 


Z(t) 


from which we can readily calculate the in vacuo modes. These modes and corresponding eigen vectors 
are given in Table 2 for data set #1. 


00 



Solution at t =0+ 


Vfe next calculate t±ie transient solution for small t , specifically for t 

For such values of t, we can ejqjloit Hie approximation: 

0 

(3/8)(^/tt) U (Ut)"^ B* 

(U/2ir)^ Cc cx)s~^c -\(l-c^))B* 


M2(t) 




Mjrt) 




Mi^(t ) 5S 0 

We note that: 


h (23L + 1) B* 

(l/4ir) [(l+2c^) cos“^c - 3c /(l-c^) J 


B* 


M3(0-h) = M3^q. 


Unlike the case for U=0, however, it should be noted that K and B are no longer zero 


< lAJ 


The Solution for large t : Steady State Soluticn 


Let us next consider* the steady state solution. For* t-^ «> , we can make the ^prxaximaticar: 


M(t) 


I - ( A + 


1 0 1 ^ 

1 

0 1 

1<M^ - )\ 

1 

(M^ - M^)"^ 1 

00. Cj 


Letting: 


A + 


(M^ - 


(M„ - M ) ^ M,(«) 

S3. 3 


1 


the asynnptotic solution at t a oo , is given by: 


(I - Ag ) Z(t) = A2 Z(t) + H u(t) + [b,Z(0)) I I 


We remark that this solution is 'exact' for U = 0. A question that naturally arises is whether tiie 
matrix (I-Ag) can be singular, in which we case we have a singular differential equation. It is 
ncaisingular for U=0, as we have alreacty seen. 



NUMERICAL SOLUTION TECHNIQUE 


For the numerical computation of the transient response due to step input, we proceed as follows. 

We take Z(0) to be zero. We choose a sufficiently small sampling interval ^ , and proceed to discretise 

as usual. Some care has to he taken because of the singularity of the kernel M^(t). 


First we note that the initial condition can be calculated from: 


M(0) y(0+) = H (as is evident by Laplace Transforming or otherwise) 


-1 


y(Q+) = M(0) ^ E 


M(0) = I - 


0 \kA 


o 


(M,-M )"^ M,(0) 
s a 3 


Next; use the discretization: 

n-1 /(k+l)A \ r _ n 

^ j J M(nA-a)da| y(kA) or + y . (kA> j ^ 


H (nA) , n ^ 1 


Let 


M 


n-k-1 


(k+l)A 

M(nA-a)do 

kA 


[In all kernels except we replace the right-hand side by the integrand at the upper limit.] 



t 


Then we have: 


n-1 


IVk-i 

0 


Since the matrix involved is "triangular", this can be solved readily without resort to iteration. 
Note that the "coefficient" of 


y(iiRLA) or is always j f^MU-a)da = M(0+). 


The transient response of the section has been calculated for the section parameters of Table 1 and 
for nondin^sional velocities of 200, 225, 250, and 275. The responses are given in Figure 12. 


00 

CO 






CHAPTER TEN 


CALCULATION OF THE GUST RESPONSE 





10.1 CALCULATION OF CIRCUIATION DUE TO GUST 

f We begin by calculating the circulation due to gust. First, we let g(x), -c» < x -< «>, 
represent the "frozen" gust velocity (field). This is a spatially homogeneous random process, 
Gaussian-distributed with spectral density p(X) v^iich can (following the Von Karman Hieory) be 
taken as: 


p(X) = Const. 


1 + 






viiere and k 2 are constants. 

From unpublished calculations, y (t,x), the circulation on the foil, is given by the 

. S 

integral equation; 


Y^(t,x) = 


- j A. 1S d? - -i r H(a,x) n (t - 0)da 

V -1 yT^ x _ e 2ir -^0 S 


vdiere fg(t) is the total circulation due to gust, defined by; 


^ Yg(t,x)dx 


to 

o> 



iO 

4=- 


and 


H(ct,x) = — nH- , z = 1 + Ua 

X - z /z - 1 

It is convenient to use the spectral representation for the gust field; 
.g(x- Ut) = /~ dG(v) 

-OO 

where 


E |d6(v)|^ = p(v)dv 


Using the fact that the order of integration in 

f :5Z1 s<f - t) a; 

-1 j^TTx •' -1 /L - ? X - ? 
can be reversed, we obtain i^on integrating in x over E-l,+lJ: 


FgCt) = 



g(C - t) d? -i- H(a) fg (t - a) da 
/I - C 2tt ■' 0 ^ 


* 



where 


H(a) = 2ir ^ -1 

A - 1 


The equation can be siraplified to; 


1 yTTuo ,, 

— I r ' (t - 0) dcT = 2 / 

2ir 0 ytJoF 6 •'-I 


yT^ 


g(C - Ut) de 


This shows that the circulation is asymptotically stationary Gaussian. We can calculate the spectral 
density in the following way. 


+2 


f dS(v) dc = (+2) r dc 


Now 


-4 •*• C g2irivc 

-^-1 


" ^-1 ZT 7 


dc - 


1 d 1 

(2iri) dv "'-1 Vl - 


dc 


iO 

cn 


= IT J«(2ttv) + ~ J, (2iiv) 
u ^ 1 


= IT £Jq(2itv) + i Jj^(2tiv)3 



vdiere we have lased: 


dx 


Jq(x) = 


J (X) = i ^ dt 

TT 0 A - 


Hence. 


= 2 f — ^ ^ g(c - Ut) dC 
■^-1 i ^ r ^' 

CO 

= (2tt) J (Jq(2itv) + i J^(21TV)) dfi(v) 


Now, substituting for H(g), we get: 

/ ^ r-(t-a)do=.S(t) 

vdiere S(IH:) is a stationary stochastic process with spectral density 


47t |jQ(27rv) + i J^(2irv)j^ p(v) 


Hence T'Ct) and T (t) are also asyrnptotically stationary, and 
S ® 


g(t) has the spectral density 


Pp (v) = 
‘g 




lL(v)| 




v^ere 


Uv) = -ir iUSdt 

J n 


2ir •'0 


yOt 


viiere 


-i- <K„(2S-) . k!(2w/U) 

2tiU ^ U ^ 


K 


/X f -St dt 


Putting X = vU, we can write 


Pp (UX) = (16w ) 
g 


jQ(2TrX) + iJ^(2TiX) 


KQ(2irX) - Kq(2ttX) ^ 


P(X) 


Ihe main result is that asymptotically (large t) the circulation is a stationary Gaussian Randcan 
Process with the above spectral density. The Spectral density depends on the assumed density for 



to 

00 


1±ie frozen gust field. We note that the latter is well apprccdjnated by the Dryden form: 


p(X) = Const. 

1 + 


especially at hi^ frequencies. We obtain this approxijiation by noting that 
1 + l^X^ , (1 + k^X^) (k^ - k^)X^ 

(T+l^xV^ (1 + k^X^)^^/^ ^ (1 + 

viiere the second term goes to zero for large X and 

1 + l^X^ ^ 

(1 + ]<2X^)^^® 1 + k^X^ 


we can also get a "timeHlcmain" representation for the gust circulation: 


I ft 

("t) = (2ir) j c^(t - a) S(a) dc 


Cj^Ct) is as defined above and has the Laplace transform! 




■1 


► . We note that F Ct) is a stationary 

O 

The circulation is cisynptotically stationary, 
c^Ct - a)S(a) da 



o 

o 


10.2 CALCULATION OF THE LIFT DUE TO GUST LOADING 


Having determined T^Ct), we can go on to calculate the lift P . Thus, we have in the 

g g 

usual way: 


Pg (-p) lu g(t) - If / ^ X YgCtjX) dx + fg' Ct) 


We need to calculate only; 


J X Yg(t,x) dx 


v\Mch, using the integral equation for yields; 

O 



X 


A - X 
■■ ■ ■■ » 

A + X 


dx 


X ~ C 



“t 1 

i-^gU - Ut)dc f f X H(a,x) dx f'(t -a) do 

1 - C 0 -^^1 ^ 


Now it can be shown that: 

rl . ,1 

J X Y„(t,x)dx = C+2) I ( 1 . - C) 

-1 g 


— g(C - Ut)d5 + r„'(t) + Ir f (Uo - v6V + 2Ucr)r'(t - a)dd 

g J 0 


Hence , 


1 "t 

P (t) = (-p)Curg(t) - 2 If / ( 1 - C^) g(C - Ut)dC - If / (Ua - 

—1 0 


It follows that Pg('t) is a Gaussian pxxxiess with the "time-<3oirain" representation; 


P (t) = (-p)EUr (t) - |r f Cq(a) S(t - a)da + — SgCt)] 
g *^g dtJ.3 dt^ 


where 


s ct) = -2 A - ^ gU - Ut)dc 

^ “V-1 


and is a stationary Gaussian process defined by 


OQ 

,(t) = -2 / 

•P- 

= -t/ 


“2iri\^t 




dc d6(v) 


e 2irivUt (j^( 27 rv) + J2(2 ttv)) dG(v) 
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Here we have used: 


1 Ji<x) 

J^(K) = J„(X) - 


= Jo(x) 


IJq(x) + J2(x)3 


and 


Jq(x) - J2(x) 


-1 -1 A - C (27ri) dv J A - C 


IT EJg(2irv) + i J^(2 ttv)] - ^ E(2ir) (Jq(2w) 


ir EJqC2ttv) - 


jQ(27rv) - J2(2m>) 


= I [J()(2iiv) + JjCaw)^ 


d? 

+ i J^(2Trv)3 


We also note that foi> large v 


J^(2ttv) + Jq(2ttv) 




1 


Note that 


/ 


tt 


- + 2t 3 dt = -i- - f ~l 


dt 


1 e® 

= -l + S_K.(s) 

s s 


= ~ ~lC(s) 

s 


= II - se K^(s) 


and hence, that 


<=3(t)dt = ^ II - 1 Kj^(s/U)3 / iK^ta/U) + Kj^Cs/U)] 
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A tine-dciiiain representation for P would be, (taking f (0) = 0) 

. S ' ■ g 

( d 

Pg(t) = (-p) j J (Uc 2 (o-) - c^(a)) S(t - a)dp + ^ S 2 (t) 


The cx^rresponding spectral density (of Idle asynptotically stationary process) is 


= ICCFourier Transform of (Uc 2 (t) - cJ(t))(l/U)(jQ(^) + (Jq(^) + J 2 (^))]|^ 


Now, the Laplace Transform of lUc 2 (t) - c^(t)J is 


= j . Kj^(sAJ) 


LKqCs/U) + K^(s/U)3 


Hence, spectral density of lift due to gust with UA = v 


111/ ^ (2ttX) \ 

= ^ K^(2irX) ) ^ i - 2i j^(2„x) ii^ pu) 


where we have used 


X (Jq(x) + J 2 (x)) = 2J^(x) 


10.3 CALCULATION OF HINGE MOMENT (M ) DUE TO GUST 

OL 


To calcailate the hinge nonant M^ due to gust, we need to calculate: 


M 


a>g 


1 12 2 
= (-p) EU j ^ X Yg<t,x)dx - aUfg(t) ~ Yg(’t»x)dx + — rg'(t)3 


Now, 


f ^ X T (t,x)dx = ^ X dx f ^ -A t , S , S , <i ,^, y? ldi; + ((1 + Ua) - Al + Uo)^ - 1) r'(t 

J-1 'g ’'•'-1 yrrif ti x-c- ■'o « 

Y ... t 

= (+2) I Al - C^)g(C - ut)dc + r„(t) - r^(0) + f c Ca) s(t - Cf)da 

J -I g g J n ^ 


li/hile 


/ ^ ■ " 2 ~' 


— (x - a)^ — - - dx f — g(C - Ut)dc + 


-1 A - c 

"t 1 

f f (x - a)^ H(a,x)dx E'(t - 0 )d 0 


Q 2ir" 


0)d0 


O 

cn 



o 

<J> 


Now it can be shown that 


Now, 


i ' f (x - a)^ H(a,x)dx 

If ~1 


( /z - 1 


1 + 2a(z - - 1) + a^(---^-- ^ - - 1) ! , z = 1 + Ua 

/r:ri I 


i /\x - a)2 . _i_ dx = i f ^ iEZI rx + (« - 2a)3dx + 

ynrir x - ? "" •^-lyr+r 


-1 


dx3 


X - c 


= -4 + (c - 2a) - (a - O' 


2 2 

= -4 - a - c + 2ac + C - 2a 


= <-P> }‘l - TgCt) - U r CO) + r_.(t) + U f 0,(0) sCt - <,)da + 2U gCc ■ 

^ ^ 0 'Ll 


-si- 

dt 


f* (X Co) r„'(t - a)da + — t t - 2ag - ; ♦ 2a) Al t ;) 

■’o ,g ^ gcc 


Ut)d 5 


-A - c 


Ut)d? 


From: 


r ' (t) = 2ir f c^Ct - a) S(a)da 
g -'o ^ 


\«e have that 


c „(a) r ‘(t - a)da = f {2a c^(a) + Cg(a)} S(t - cj)da + (i + a^) S(t) - * 

J 0 12 g *^0 

since (as can be verified by taking Laplace Transforms) : 

(i + a^) - c,(t - o)da = (i + a^) 6(t) 


Hence , 

2 "t 

M (t) = (-p)((l - a)U r (t) - U r (0) + ^ - - r_'(t) + u f CgCa) set - a)da -.U S^Ct) 

a,g. ' g g ^ g 0 

^ (2a c^(o) + c^(a)) S(t - a)da + r'(t) + — J - 2ac - C + 2a) ^ gU 

Ql ^ q ■ d"t "“X ^X Q 


Ut)dC> 


O 
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New 


f (C^ - 2aC - ? + 2a) - ----- d? = [ (2a) A - e2TriCV _ f ^ kA~^ ^2TrxCV 

-'-1 yT^ ■'-1 ■'-i 


= TTa rJr,(2irv)) Hh J_(2irv)] — (^) ~lJ„(2irv) + J^(2irv)] 

^ ^ 2TTi ^ dv ^ / 

J„(2ttv) 2JT(2Trv) J«(2Trv) 

= TTa lJ„(27rv) + J2(2irv)J - iri = (ira) - iri — 

2itv 2tfv 2ttv 

where we have used: 


J«(x) 

j'(x) = -2 -H J^(x) 


J*(x) = -J^(x) 


2J-, (x) 

4— = J„(X) + J^(K) 


This shows that ^Ct) is asynptotically a stationary Gaussian process (if we may set TgCO) = 0) with 
finite second ironent. The spectral density can be deduced readily fratn the timerdomain '‘convolution" 
integrals. 



10.4 CALCULATION OF FLAP MOMEOT DUE TO GUST 


We have 

m| = (-p) I f U(x - c) Y„('t>x)dx + — f (x - c)dx f y Ct,y)dy | 

^ r c S dt ■'-1 2 ) 

= (-P) jf^ U(x - c) Y„(t,x)dx + r^'Ct) - ^ S t -Sl i y (t,x)dx 

C-' c ® 2 ^ dt ■' c 2 ® 


f (x - o) Y.(t,x)dx = (-1) f (.h (x - o) — ~ ^ dx f - i - * i SiS-T-ffil 
■'c e > ir ■’.l.TTc x.c 

1 f ^ r ^ — V 

"T I Cx - c) ^ H(a,x)dx f '(t - a)ckj = + T„ 

ir^ J 0 I c S 1 2 

The only way to handle the first term T^ would appear to be to integrate with respect to x first, and 
then use the spectral representation for the gust . The second term T 2 , of course , presents no problon. 
Thus, we have for the integral with respect to x in the first term: 


(-1) (|) (x - c) • — i— dx=|r^[x+(C-l-c)+ — 


/I + X X - C 

r — I X - + (C - 1 *- c) OOS~^C + — ^ ^ ~ ^ ~ • log 


1 
c 

2 

A - 


C ~ C - cc 
X - C 


dx 


1 -x2 


/a - cKl 7~cT t /(I + c)(l ~ C) 


/(I ^ c)(l + c) - /(I + c)(l - C) 


o 

lD 
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viiere 


g^(v) dG(v) 


g^(v) = I / ^ rc/r. o' 


(1 + c) oos'^o + _£1 log I . |3 !^ ., t ., <J5 

/l - c2 A - c 




Cl + c) cos ^c) + cos ^c3 — — — • " ■ dC 

AT::T 

+ “ f ^ g2TTivc £(c ^ 5) log |.p 


Only the second integral requires attention. We use integration by parts to sinpli:fy it. Thus 


a_ log 1 1 - • 1 

K 




r c - c 


/: 


(o 


- e)dc = 


c e 


2iviv? 2Tf±xr, 


^2irivc ^ 


27Tiv (27Tiv)2 C2iriv) 


g2TTivc 

(2Tfiv) 


Kg - c) + — ^ ] 

2'iriv 
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Next in T^: 


.t 


1 f f * 

-J I I - c) - ^ H(tf,x)dx r '(t - a)dcT 

7T 0 C yi + X ^ 




V ' 0 


vheve 


H,(a) = (.4) (z - c) tan'^ Az. + PCI - c)_ _ ^ + 1 

-L /7— — : r- =- 


- - 2 [(1 - z + c) cos’-^c 

r Vz - 1 


Next we need to calculate: 


d f ^ (x - c)^ d „ d 

— I Yp(t,x)dx = T + H„ T 

dt c 2 e dt ^ dt ^ 


T = r2\ f ^ (x - c)^ A ~ X 
^ ’ \7 C 2 A~n^ 



Now 


Hence 


J (x - c) 


2 A - X . _J 
A + X X - ? 


dx 


= f C-x^ + (1 + 2c + C)x + (C + 2c + - c^ - 2c)3 a + f 

^ c A ~ X- ■' c 


(c ^ c)^g - 

>4 - x^ 


= HU (c) + ~ log 

^ A - 


Al - cKl + C) + Al + cKl ^ C) 


Al ^ c)Cl + C) - Al + c)(l - C) 


Now let 


Tg = - ^ J m^(c) g(C Ut)dC " 7 " Ut) log | j dc 


f\c - log I 1 dc = fS; - o)^ log j j d? 

•^-1 -1 


27rivc 

e 


(c - 


-1 


„2irivc 

e 

d2 


(2iri)^ 

CM 



/ log I I dc 


C) dx 
X - C 


CD 
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Hence it is enough to calculate 


/ ^ log I I d5 = 

-1 2'rriv 


Hence, finally. 




^~2irivUt 

e 


EhgCv)^ dG(v) 


where 


-i /r^ 1 C2wi)2 dv^ 


^•^2'irivc f(c,v) 
2iriv 


Nejct 


^'^4 " Tg’Ct - a)da 


dt 


vdiere 



Hence, finally. 


g(C - Ut)dC - ~ f r_'(t - 0)dcf 

2tt 0 ® 

og I 1 g(5 ■*' Ut)dC 


+ L.±.f^ H (a) T »(t - a)da 
4 dt-^0 2 S 

vftiich is tfie time*-doinain version of the flap nonent due to gust. It is seen to be Gaussian, and asynp- 
totically stationary, vdth spectral density that can be readily deduced fran; 




cn 
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M|(t) 


= (-P)< 


00 

-.00 


f(cjv) + 2(»4 - - (1 + 


+ 2(cos’'^c) 


\2Triv, 


Jq(2ttv) - 


c)cos''^c)(jQ(27rv) + i J^(2wv)) 


CJ^(2irv) - 2iriv J^(2 ttv)) 


2ttv 


^-2TrxVUt 


- — f r '(t . 

2 u ' 0 ^ 




g 


f HAa) r^At ^ a)dcr + i [ h.(v) dG(v) 

HIT dt 0 • S ^ -Loo ^ 


vtiere 


^2irivc ,2 

hsCv) = 2 2- S' 


(2Tri)"^ d 


^-2'irivc f (c,v) 
2'iriv 




2irivc + X 

/l^ 


(C + 2cC + - 2c - i)cos''^c + / - (1 + ~ + 

2 


a)dCT 


C) dc 


10.5 CALGULATICai OF STRUCTURAL RESPONSE TO GUST 


The structural response is governed by the dynamic equations; 


M X + B X t K X = 
s 8 B \ jn 



P„(t) 

g 

M^Ct) 

M|(t) 


where P^Ct), M^(t), l^(t) have already been calculated. Prom these calculations > we can see tiiat 
the right-hand side can be e^jpressed as : 

/“ Uv) dG(v) 

*^-00 

vAiere, however, the transfer function L(v) is not necessarily "physically realizable". In other 
words, we cannot egress it as; 


L(V) = /' 


-2irivt 


W(t)dt 


where W(*) would be the *V?eighting matrix" corresponding to a physical system (vanishing for 
negative t— values). The gust itself can, if we follow the Bryden form of the spoctral density, 
be expressed as: 


-<1 
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A P P E N PI X ONE 


INTEGRAL FORMULAS 
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Forroulas 




'-1 


dC 


TT, 


5-1 •;£:? 


x-c 


dc =0 |xj < 1 (cf.ba) 


2 . 


.O+C dC 


-1 


x-«. 




/C? 


x-«. 


(-1) rn (x+1) 


'-1 




Cl - 3dc = - IT, |xl < 1 


3 . 


^ X, 


-1 


x-C 


dC = 


f V XT 

J-1 /r:: 


l+C r 


C— -p- -l]dC = = irx - X 


u. 


_L 

5-1 


vT?c 1 




C-z 


dC 


Nov; C -i~ dC 


J 

1 


1 . 
/ T 


[i^X]dc 

•■C-z ^ 


cos <J» = C 


^ Cl + ^3dc 


vil? 


4> + (z+l) 


G 


A-c^ 


JLl 

C-z 


dC 
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/i-C' 


1 ,, 
C— z 


1-t^ 

l+t*^ 


i-c' 


= 1 - 


(1-t^)^ 

o»?? 


4t 


(i;-t2)2 




2t 


(1+t^) 


hi 

i+c 


dc 


((in-^)(-2t) - (i-t^) ( 2t ) )'d t 


(1-^tv 


_ - 4t dt 
(1+t^)^ 


..1 -. 

5 ■— dc = f 


2t 


(l-t^)-z(l+t^) 


dt 


f. 3-+C 

(-1) J (2) 


0 


dt 


(z-i)+r(i+z), 


If z > 1 


=> 


’'17^ 
(-2) J 


dt 


0 t^(l+z)+(z-l) 


Cn 


V(l-c)/l+c 


l+z 


'I \ 

t y 


/2-^i Tan 

Vz-£ 


-1 


‘'Z'i 

/z-i 



1 




(-2) » /n^)>CT' 

^“2” /O+c).^ 


rf 2 < 1 , 



(- 2 ) 

1+2 





/TTc 

rr:? 


^0 


Log 


t~A-z//l+z 

t+vA^/ZE+z 



dC 

C-2 


1 /n:^»A+2 +/TI+c)Ali‘ 

^_^2 /(1-c) -/(l+c)A-z 



CO 


dt 

XI^ 

1+z 


, |z| < 1 
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6(a.\ 


1 

5. 


' 2 


/i- 


~ d.; 


- Q 2 < 1 


Hence 


1 

7 C 

/I+t 

1 

3 

c 

»4-£; 

i-z 

8. j" 

•/DT 


c 

AT 

C-z 

1 

9. C 

A+c 

c 

^-1 

AT 

C-z 

1 

10. \ 

A+c 

1 

Al 

/!< 

C-z 

10a. ( 

AT 

(C-a) 

J-1 

AT 

x-C 

10b. 1 

/IT 

(C-c) 

J-1 

AT 

x-C 


dc = 


_-l ^ /L+z , 

^ /iZi* Log 


Al-c) A+z -»^(l+c)A-z 


, z‘ < 1 


-1 


^z+i' Ai-c) A+l 


, z > 1 


dC = IT ; Z < 1 


= IT - IT •, Z > 1 


-1 


r^C + log I 1 .1 (x-c) 

vToc 


J 1 



LL. 


‘ -1 


__jl_ 

C-l-Uo' 




x-C 


✓ 17c 

71^: 


^■ 1 . 1 . 

y.-z *c-3-Uj x -c 


ir 

y-7 


T1 

x-z 


JL 

x-z 


/z»^l 


IT /z+1 


12. H(o) = - C 

TF J 


7T 


[|] 

U 


S‘ 

/1=^ 

H(a,x)dx 


(-Tf) 

✓z+r 


/T=^ 


/z-I 

'-1 


(-TT) 

.ct: 

1 

‘ $-1 

1 

.CT 

/!-✓ 

vOT 

^(l-z) 

V?: 

(-2) IT 

- ■ 2 

1 

- tt) 


X-Z 


dx 




to 

cn 


(2ir) (1 + 






(27f) (,/ZTT -1) 


✓?=r 


= 1 + Ua 
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13. 


C 1 

J_1 /2+^ 


-1 


(I4 

x-C 


- 1 )clx = 


“ C < 1 


>1-X 


13(a). 


C ^ 1 

'-1 ^**2 ^ = 






/'l“Z 

^rrr- - Ddx 

/V— ^ 


J(z-l) 




- 7T, Z > 1 


I'l. 


-L 


✓17c 

/R 


Cdc 


tt/2 

V (l+sin0)siji0d0 

-TT/2 


TT 

2 


15. J ^ ac 

o A'Tf 


.-n/2 

j (l+sin9)d0 


“ ^2 ~ ■*■ ✓l-c^ 


= cos c + 


• ✓TT? 


r ^ c dc 

C /Itf 


TT/2 

O>uin0)si,OuO = .1 


r ^ <^-c)dc = . cl - 


c ./ir 


c) go; 


^-1 

3 c = 


= + (i-c)cos“\ 


t is- 


1 i 
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J6. 


f 2-z; 


/L+C 2, 

C dC 


n/2 


J (lK*;i.Ln0)3in‘^0dG 




sin.j* = c 


tx/2 


I 


.2, 


ir/2 


sin OdO + f sin^BdG 


- J cos~^c + i c/C? + 


c 


Ti/2 ^ 
C OS *^9 


17. 


“1 A+x 


X dx = 


f (l-x)x J 9 

) , *7=^ ^ f sin^ede 

-1 a :7 ^-”/2 


TT 

2 


18. J (-i+2c-2c^)dc 

•'c ATc 


ir/2 


= C (l+sin0)(-l+2sin0-2sin^0)dO 
•’‘I* 


-J^/2 

= 5 (-l+2sinO-2sin^0-sin0+2sin^0-2sin^0)dO 

(I* 


COS0] 
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j (•'i+s.lriO- 2sin^O) JO 


= -cos 




(-l+2c-2c^)i;;dc 


c /l-c 


^ (-sin0+sin^0-2sin^0)d0 

fh 


(-1 ) /L-c^ + i cos 


I cos“^c + ~ 0/C7 


("1 cos-V 


18 (a) 5 ^ m (-1*25-252 ^ _ 

-1 i4-c 


19(a) J ^ (-1*25-2 c')W5 = - 1 
-1 /l^ 4 


3sin20 _ sin^0cosO-, 

16 “ u 2 



20 , 




( a - c)dx “ 


if/2 


(1-sinO) (r>i;i6-c)d0 


r 

= (1+c) A-e 


-1 i -1 1 /T2 

C cos c - 7 j- cos C - 7J- C/l-G 

if ^ 


= (D.+ - (c + i ) eos ”^ G . 


21 . 


-1 yrHv « q 

\ x(x-c)dx = \ (( 1+c) sin 0- csinO -sin'^0)dO 

c A+x 


,1+Cx (l-cos^0) 




- csin0-sin'^0)d0 


4> 


1+c 

[ ( i ~)0 - 

<!> 


3 

(1+c) • 2» . Q cos 0 . 

-q— sin 0 + CCOS0 g — + COS0J 


1 — 

Y (1+c) COS "^G - -~2 


. .,, 4:7 


^ (1+c) cos’*^c + y A-c*^ (c^-c-2) + 


(^ 


2^3 


N 3 

to 



GJ 

O 


\ ~ (X-l)(X“C)cjx = 

•'c /I+x- 


c -j y — /• / _2v3 

- ^ cos c + ^ (.4-c") + ^LirSjL. 

^ 3 


23. 


7T/2 


S (x-c)^dx = C (s.-mQ-c)^ (1-Sin8)d0 

_i 


c /L+^ 


tt/2 


[ ^ (c2 t 1 - c)6 - IMhI 


sln2a t (l+o)^ oose - s°5_® ] 


= c| (U2cnc^) cos-^c + 


= [ 


(l+2c+2c^) 


cos"^c + ~ (-2c^-9c-4)/C?] 


24. 


= h^(c) 


/T^x 2 

3 g — (x+l)(x-c) dx = C (sin0-c)^(-sin^O + l0 

c /1+x 


= h^Cc) 


25 . 


-1- "H^/2 

( ~~ dx = f (l-sinO)cIO = cos"-^G - -/C? 

c /1+X' •' ^ 


26. 




2^ J_ 

c 


c3x 


i' 


,/[:7 


C -1 + ^]Jdx 

X-2 


= cos 


-Ic 4 Tan-1 iks.).^ . , > 1. 

(l+c)/z-l 


/z^] 


27. 


I. ^ 


/U5T 


‘li/2 ■■ X - 

dx = f (sinO-sin^O)dO = /Z? - i cos“^c - 

2 2 


c /l+x 


28. 

f, (t) = ^ 
1 If 

(x-c) dx 1 



-Ic ✓I+x Jc 2llc 

29. 

f.(c) = i 
2 IT 

f" (x-o)dx dr 

2i:+3r Jc 


^ X 2 —1 

'C cos c 


(c^- 


-»■=(- i) { (cos“^c + /C?)[(l+c)cos"^c + 
- 2)/3?+ 2 (/I^)33 ^ 1 


c 
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al|M =9|H =S1H- 



o 

CO 




APPENDIX TWO 


SPECIAL FUNCTIONS 






where we have used: 


)o a::? 


d? 



Now the second integral in brackets can be integrated by parts noting 
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Now the first integral; 


I 


A-x 
c 


(x-c)dx 


= 2 C(2+c)v^-c^ - (2c+l) cos’^cJ 


Hence 


= i C(2cHKoos*^c)^ - (2 +c)/{? 03 s- 1 c - + cAJ 


COS~^c] 


I C(2c+1Kgos“^g)^ - (1-c^) „ 2/C?cos”^, 


cos ci 


Next 


f,(=) =| f" j 52 

*'c /I+x ’c A!^ 


. ^ 


d? 




TT 


•2' f" iS 

■'c /L+x Jc 

= (2£, f2- . 

). J %id; 


/ ^ a- d, 


✓TR* 


c /T+x 


c 


* 4) f" .„2 


Cx'"-Gx)dx f 
-'c 


dS 

X-? 


The, first term is 


= (^) CCl+c)cos~^c + (/C?)(c^-c-2) + i . 


) 3 • (cos”^c + 
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The second term is: 



✓TO 



-cx) 


✓1^ 


log I 


dx 


and by integration by parts 




J A~c^ (2 oos^e - c cos6 - c^]d9 


since a ready integration shows tha.t 



Hence 2nd term 

■ “ i? Ccos"^c + c/C? - cX? - c^ cos~^c] 

= “ ^ Jl-<^ (1-c^) COS~~C 


Hence 


faCc) = {C(l+c)cos“^c + (c^-c-2)>/C? + i(/C?)^i(cos“^( 


c+ 


g^(c) 


. 2 
tf 


(x-c)^ ^52 dx 

Vl?? 1 

II 

'c 

✓1+x 

vll? 

. ^ 

IT 


(x-o)^ dx [- 

oos c - — — 

II 

^ c 


/1-x 

. 1 

^ IT 

(- cos“^c) f (x-c)^ 

— dx 

II 





2 

TT 

j (x-c)^ log 1 

1 dx 




log j (3 dx 


2nd tera, by intagration by parts 
. Cx~c)^ 


1 

x-c 


dx 


ax 


31T j, 


■la’ r Ccos^S ■<• - 2c GOs0)dS 


* c/il? . c W^c - 2=^] 
2 


d7T 


= - —• C(l+2c^) cos”^c-|<=^^? 
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The first terra 


2 

ir 


2 

if 


Hence 


g^Cc) 




g2(c) 


COS c 


-1- (x-c)2 ^ dx 


cos'^c [( ^r i yj £, ) cos“^c - i (2c^ + 9c + 4) /I? ] 


. / >^T^ rCy'Cj (l+2c^) ^ 1 ,^_2 ^ ^ 

" ^ 3 ^ cos C + (2c + 9c + 


4)] 


I (cos"^c)^ (l-2c+2c^) 


(-2) 



(x-c)^ dx 

^'c 

/TfiT 

■ 1 



Ic ^ 


/ v2 /1-x 
(x-c) < 

^ c 

vnj^ 


(x-c, 2 ^ X 

^c 



i: 






X-? 


d? 


dx C 


iEi n > 




(1 > ^)do 


/I=Q 


dC 



Now 


\ m: 1 


cos c + — — log 

✓1-x 


Hence 2nd term in g 2 (c) is: 

- I y x(x-c)2 ^ dx oos"^c 

yl+x 


2 5 

~ I xCx-c) log 1 1 dx 
•‘c 


s - h,(c) - I 


.X* cV 

c^+ <= X 
C 


f 


2cx^ 


) iS 


dx 

x-c 


Now 


hjCc) = - 


,2 


f xCx-c)“ dx 

/I+x 


- ^ ( (cose - c)^ cose (1 - cose) de 

'0 

2 r ^ 4 

* 7 in ® + 2c) cos^e - c^cose) d6 
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- 2 ; 3 ^ 3 . 1 ' O .3 

"7^78“^ ’a “ 7 ^ * (l+2c) (sin.^ - ^ 

0 

- (c + 2c) ( I* + j sind oos((>) - sincj)} 



= i {(ccs-^c) (o + - |) * - 5^ 


c)} 


Next 


" ,=y 


i f (i- 


2cx o 




w 


dx 

x-c 


= I (~) (3x^ - 5cx^ + c^x + c^) 


3n dx 




■fe-c? C 


(3 cos^9 - 5c cos^9 + c^oos8 + c^)d9 


= i- { - I- cy^? oos”^c + 2 + i ^ 


6ir 


5 2” 


Hence 

g^(c) = (cos ^c + 1 ^-c^) (“ C(l-2c + 2c^)cos ^c - y (2c^ + 9c + 4 ) Jl-<^ 

- I C(c + |i- f) cos"^c + - 

- ~ Cc^ _ y C OOs"^C + (2 + - ^ )] . 
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TABLES 


<> 




TABLE 1 


Section Parameters 


= 0.2 




= 0.0125 




= 0.25 


0) 

=100 rad/ sec 


Ct 


= 0.00625 



= 300 rad/ sec 

= - 0.4 



= 50 rad/sec 

= 0.6 




= 1 




m 

m 



« S' 

s 

= 40 


Trpb^ 

TTp 


P - 1 





m ~ (4 q)tt 
s 



0 


TABLE 2 

System matrices ; U = 0 ; ? o 

p 



O.lOOOOOE+01 

0.200000E+00 

0.125000E-01 


0.200000E+00 

0.2S0000E+00 

0.187500E-01 


0,125000E-01 

0.187500E-01 

0.625000E-02 


K 


0.250000E+04 

O.Q 

0.0 


0.0 

0.250000E+04 

0.0 


0.0 

0,0 

0.562500E+03 





-0.250000E-01 

-0.999999E-02 

-0.580566E-03 


-0.999999E-02 

-0.712499E-02 

-0.473441E-03 


-0.580566E-03 

-0.473441E-03 

-0.S57130E-04 


0.0 

0.0 

0.0 

^3x3 

0.0 

0,0 

0.0 


0.0 

0.0 

0.0 


■0.293400E+04 

0.251430E+04 

-0.3SS226E+03 


0.251430E+04 

-0,147478E+0S 

0,894246E+04 

^3x3 

0.157874E+04 

0.397441E+0S 

-0.115730E+06 
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TABLE 2 : ( continued ) 



0 0 

^3x3 -.01485 

-.0006165 -.0005549 


0 

-.00245 

-.0000917 


(I - A3)“^A = 



-.002943 2465.6 -451.0 

.002457 -14636 • 3676.0 

-.05415 45319 -104830 


1 3x3 
^3x3 


TABLE 3 


I"fodes and Eigen Vectors for U 
Eigen-values + 330.22 j 

Eigen-vectors : 


1.9 

X 10“® 

-2.8 

X 10“^ 

3.01 

X 10”^ 

+ 6.33 

X 10“^j 

+ 9.24 

X 10"^j 

+ 9.257 

X 10"^j 

Eigen-values + .106 . 2 j 
Eigen-vector 

+ 2.357 

X 10"^ 

- 8.75 

X 10“^ 

- 4.25 

X 10“^ 

+ 2.50 

X lO'^j 

+ 9.297 

O 

I 

+ 4.51 

X IQ-^j 


Eigen-values + 47.985 j 
Eigen-vectors 



2.00 

X 10’^ 


5.64 

X 10“^ 


2.39 

X 10~^ 

z 

9.61 

X lO’^j 

z 

2.71 

X 10“^j 

z 

1.15 

X 10”^j 


FIGURES 














































































































































































1. Report No. 2. Government Accession No. 3. R^ipienft Catalog No. 

NASA TM-81351 


4. Title and Subtitle 

CALCULATION OF THE TRANSIENT MOTION OP ELASTIC AIRFOILS 

6. Report Dan 
August 1980 

FORCED BY CONTROL SURFACE MOTION AND GUSTS 

6. Performing Organization Code 

7. Author(s) 

8. Performing Organization Report No. 

A, V. Balakrishnan and John W, Edwards 

H-1125 


10 Work Unit No 

9. Performing Organization Name and Address 

505-36-24 

Dryden Flight Research Center 


P, 0. Box 273 

11. Contract or Grant No. 

Edwards, Calif, 93523 

• ' 

NAS4-2472 


13. Type of Report and Period Covered 

12. Sponsoring Agency Name and Address 

Technical Memorandum 

National Aeronautics and Space Administration 


Washington, D, C. 20546 

14. Sponsoring Agency <3ode 


15. Suppiementary Notes 


This work was performed partially under NASA contract NAS4-2472, 

16. Abstract 


The time-domain equations of motion of elastic airfoil sections forced 
by control surface motions and gusts are developed for the case of incom- 
pressible flowc Extensive use is made of special functions related to the 
inverse transform of Theodorsen's function. Approximations are given for 
the special cases of 25 ero stream velocity, small time, and large time, A 
numerical solution technique is given for the solution of the general case 
and examples of the exact transient response of an airfoil are presented. 


17. Key Words (Suggested by Author($M 18. Distribution Statement 

Flutter 


Active control 
Aeroelasticity 


Unclassified - Unlimited 


Subject category: 02 


*For sale by the National Technical Information Service ^ Springfield^ VA 22161 


■>r 



19. Security Classif. (of this report) 

Unclassified 


20. Security Classif. (of this page) 
Unclassified 
















